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Theoretical fit of Cepheid light an radial velocity curves in 
the Large Magellanic Cloud cluster NGC 1866. 
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ABSTRACT 

We present a theoretical investigation of multifilter (U,B,V, I and K) light and radial 
velocity curves of five Classical Cepheids in NGC 1866, a young massive cluster of the 
Large Magellanic Cloud. The best fit models accounting for the luminosity and radial 
velocity variations of the five selected variables, four pulsating in the fundamental 
mode and one in the first overtone, provide direct estimates of their intrinsic stellar 
parameters and individual distances. The resulting stellar properties indicate a slightly 
brighter Mass-Luminosity relation than the canonical one, possibly due to mild over- 
shooting and/or mass loss. As for the inferred distances, the individual values are 
consistent within the uncertainties. Moreover, their weighted mean value corresponds 
to a distance modulus of 18.56±0.03 (stat) ±0.1 (syst) mag, in agreement with several 
independent results in the literature. 

Key words: star clusters - NGC 1866 - Cepheids - Variable stars . 



1 INTRODUCTION 

Classical Cepheids are considered the most important pri- 
mary distance indicators within the Local Group. Their 
Perio d-Luminosity relation, dis covered by Miss Leavitt in 
1912 l|Leavitt fc Pickering! [iQil ') for Cepheids in the Small 
Magellanic Cloud and usually c alibrated in the Large Mag- 
ellanic Cloud (LM C) (see e.g. iMadore fc FreedmanI [l99ll : 
lUdalski et alii 19991 ). is n ow at the basis of an extragalac- 
tic di stance scale (see e.g. iFreedman et al. l l200ll : ISaha et al.l 
I2001I . and references therein). Indeed, with the capabilities 
of the Hubble Space Telescope, Cepheids have been observed 
at distances (up to ~ 30 Mpc), enabling the calibration of 
several secondary distance indicators capable to reach cos- 
mological distances and t o provide an estimate of the Hub- 
ble constant Ho (see e. g. IFreedman et al.ll200ll . for detailed 
discussion). In spite of the most recent relevant efforts in 
the direction of reducing the uncertainty on the Cepheid 
based extragalactic distance scale (see e.g. iRiess et al.ll201ll . 
|2012| . and references therein), some systematic effects, in- 
cluding the effect of the host galaxy metal content, re- 
main unsolved, and different authors keep to provide signif- 
icantly different estimates of the Hubble co nstant (see e.g. 
iTammann fc Reindll |20]1 : iRiess et al.ll2012l . and references 
therein). The first crucial step for the calibration of the ex- 
tragalactic distance scale is the distance to the LMC. Several 
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methods have been adopted in the literature (see e.g. lWalkerl 
I2OIII : iMolinaro et al.l |2012| . and references therein) provid- 
ing values that range from ~ 18.1 to ~ 18.9 mag. Systematic 
effects such as a non negligible metallicity dispersion, differ- 
ential reddening and a significant depth of the Cloud, are 
known to be at work. 

In this context, classical Cepheids belonging to young 
stellar clusters in the LMC play an important role, being 
at the same distance and sharing the same age and chemi- 
cal composition. Thanks to these advantages they offer a 
unique opportunity to investigate the uncertainties affecting 
both empirical approaches and theoretical scenarios. NGC 
1866 is one of the most massive young clusters in the age 
range 100-200 Myr, and it has been the subject of a very 
long list of paper s, star ting w ith the pio n eering ones by 
lArp fc Thackeravl (|l967l ) and iRobertsonI l| 19741 ). Subse- 



testbed of stellar evolution theory (e.g. iBrocato et al. 


1989, 


1994 2001 12004; iBarmina et al.l 20021: Walker et al. 


2001: 


Chiosi et al.1 Il989l: iTesta et al. 19991). as a Cepheid host 


fWelch et al' 


I991I: Gieren et al.' 'l994': 'Welch & Stetson 


1993: Walker 


1995: Gieren et al. 2000: , Storm et al. 


2005|: 


Testa et al. 2007: Molinaro et al.ll2012h. or as a dynamical 



of a strong debate over the presence of convective overshoot- 
ing in intermediate-age stellar models, and on the fraction 
of binaries in the main sequences. Moreover, this cluster lies 
in the outskirts of the LMC, so that field contamination is 
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not severe. The investigation of Cepheid properties in this 
cluster can provide crucial information for our understand- 
ing of the physics and the evolution of intermediate mass 
stars. In particular the comparison with the predictions of 
pulsation properties based on hydrodynamical models is an 
important tool to constrain the individual distances and 
the intrinsic stellar parameters, without relying on stellar 
evolution models, and in turn independently constraining 
their physical and numerical assumptions. In particular, 
the possibility to obtain this information from the direct 
comparison between mod eled and observed light curves 
has been first cla imed bv | Wpod et al.l l|l997h for Classical 
Cepheids and bv iBono et al. I 2000ah for RR Lyrae. Sub- 
sequent applications both to field and cluster pulsating 
stars, sometimes including the additional match of radial 
velocity or radius curves, have provided self-consistent 
results, also in a greement with independent estimates in 
the literature (seelDi Fabrizio et"al]|2002l: iBono et al.ll2002l: 



Marconi & Clemcntini^ 
2007 ; N_atalo ct al. 200i 



arconi Degl'Innocentil 
iMarcorii ct al. 2oTol v 



12005 



In this paper we present an accurate comparison be- 
tween observed multifilter light and radial velocity curves for 
a sample of Cepheids in NGC1866 and the theoretical coun- 
terparts based on nonlinear convective models. The photo- 
metric and radial velocity data are introduced in Sec [2] while 
the fitting procedure is described in Sec[3] The results of our 
analysis are contained in Sec[3]and include the best fit struc- 
tural parameters, the comparison with spectroscopic data, 
the implications for the distance, including a critical discus- 
sion of the associated uncertainty, and the Mass-Luminosity 
relations. Finally, Sec [5] contains the conclusions of the pa- 
per. 



2 THE DATA 

The adopted photo metric data include observations in the 
U, B, V, I bands (|Musella et all |2006| . Musella et al. in 
preparation) with t he addition of the near infrared K band 
l|Storm et al.ll2005l : Westa. et al.ll2007h . To properly sample 
the region near the maximum of light, we have integrated 
the B, V, I photometry of the Ceph eids HV 12198 a nd HV 
12199 usin g the observat i ons fro m .Welch et all (ligQll ). MA- 
CHCQ and |g ieren et al. I (|2000h . Moreover, in some cases 
we excluded one photometric band because the light curve 
was poorly sampled and/or lacking the maximum and/or 
the minimum phases. 

To compare the Near Infrared data with models we have 
transformed the K band measurements from CIT and LCO 
into the Johnson photometric system, using the relations by 
iBessell fc BrettI ([l988): 



Kj = KciT/LCO + 0.027 - 0.007(V 
As for the radial velocities. 



CIT/LCO 



(1) 



w e used data from 

IStorm et al.1 ll2005h IStorin et all (|2004 ). IWelch et al.l (Il99ll ) 
and Molinaro et al. I l|2012l ). 

Table[T]summarizes the adopted number of photometric 
and radial velocity measurements for the selected Cepheids. 
We phased the light curves by requiring that the B band 
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Table 1. The period P, the apparent V magnitude and the num- 
ber of measures for the U, B, V, I, K hght curves and radial 
velocity curves of all selected Cepheids. 



Name 



Period V{mag) U B V I K Rad. Vel. 



HV 12197 3.143742 15.91 

HV 12198 3.522805 15.77 

HV 12199 2.639181 16.09 

We 2 3.054847 15.86 

V 6 1.944252 15.97 



3 69 87 38 35 38 

30 69 90 62 77 38 

34 157 199 172 54 39 

15 69 90 62 5 12 

21 69 90 62 10 10 



maximum of light occurred at phase zero, while the maxi- 
mum in the other bands were shifted in phase as expected 
(see e.g. iLabhardt. Sandage fc Tammanr] 1 19971 : iFreedmanI 
Il988l ). Similarly, the radial velocity curves were phased by 
requiring that their minimum occurred at phase zero. 



3 MODEL FITTING 

New nonlinear convective pulsation models have been com- 
puted to reproduce the observed luminosity and radial veloc- 
ity variations. The adopted theoretical framework is based 
on a nonlinear radial pulsation code, including the non- 
local and time-dep e ndent treatment of turbulent convection 
JStellingwerj Il982l : iBono fc Stellingwerj 1 1994 iBono et al.1 
Il999l ). The system of nonlinear equations is closed us- 
ing a free parameter, a^j, that is proportional to the 
mixing-length parameter. Changes in the mixing length 
parameter affect, as expected, both the limit cycle sta- 
bility (pulsatio n amplitudes) and the topology of the in- 
stabi lity strip l|Fiorentino et all l2007l : iDi Criscienzo et al.l 
|2004 ). Similar approaches for the tr eatment of convec - 
tive transport ha v e bee n developed bv iFeuchtingerl l|l999l ): 
iBuchler fc Szabd (|2007l ): lOlivier fc WoodI l|2005h . We re- 



mind here that the nonlinearity and the inclusion of a non lo- 
cal, time-dependent treatment of convection and of its cou- 
pling with pulsation allows us to accurately predict all the 
relevant observables of stellar pulsation, namely the com- 
plete topology of the instability strip for each selected pulsa- 
tion mode, the accurate morphology of light, radial velocity 
and radius variations and the associated pulsation ampli- 
tudes. 

The modeling of the observed light and radial velocity 
curves has been organized in three main steps; 

• First, we constructed a set of models with fixed 
chemical composition (Z=0.008, Y=0.25, consistent with 
the abundances meas ured for Cepheids in NGC 1866 
(|MucciareUi et al.|[201lh ). mass and period (equal to the ob- 
served one), varying the effective temperature and luminos- 
ity in order to reproduce the observed variations. 

• Once identified the best effective temperature from the 
previous step, we built a sequence of models at fixed chemi- 
cal composition, period and effective temperature, by vary- 
ing the mass and the luminosity in order to obtain the best 
fit model, reproducing simultaneously the multifilter light 
curves and the radial velocity one. 

• In some cases we also changed slightly the metal and 
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helium abundance within current uncertainties on the LMC 
chemical composition. 

Each model curve was phased in order to find the max- 
imum of light in the B band at phase zero. Afterward, 
for each model we calculated the shifts in magnitude and 
phase, which gave the best match between the theoretical 
light curves and the observational data. Specifically, for each 
model, we calculated the phase shift, Scj), and the magnitude 
shift, 5M, which minimized the following function: 



2 

X = 



E E 



+ S<i,') + SM')Y (2) 



where the index i runs over the number of photometric 
bands, Ntand, and j over the number of measurements, 
Npoints- In the equation above we used a spline interpo- 
lation to evaluate the theoretical absolute magnitude Mmod 
at the phase (f)j, of the jth photometric measurement rrij, 
plus the shift 5(f). We note that the parameter 5M^ gives the 
distance modulus of the analyzed Cepheid in the ith photo- 
metric band. 

As for the radial velocity curves, we transformed the 



model pulsational velocity. 



puis 



into radial velocity, u' " 



by using the projection factor p in the relation u''" — 
_ij^p"'s and phased it in order to find its minimum value 
at phase zero. In the previous formula we fixed the pro- 
jection factor to p=1.27, obtained by using the mean pe- 
riod for the selec ted Cepheids (first-overto ne fundamental- 
ized according to iFeast &: Catchpoll 19971) in the equation 
p = 1.31 - 0.08 logP. bv lNardetto et al.l (|2009l ). 

Then, to match the model radial velocity variations 
with the data, we minimized the function of an equa- 
tion similar to eq.Q, where the magnitude shift is replaced 
with the 7 bary centric velocity, namely: 



RESULTS 



puis 



+ ( 



+ 7 



(3) 



In the present section we discuss the results obtained from 
the fitting procedure. The match between the best fit models 
and the data is described for all the selected Cepheids. The 
derived structural parameters for all Cepheids are given and 
for some of them we performed a comparison with the results 
obtained from the spectroscopy. Finally, we will give the 
best fit distance moduli and reddening values and discuss 
the implications of the derived Mass-Luminosity relation. 



dex in log L/Lq). Two of these "secondary" models have 
the same temperature of the best fit one and a varied value 
of the luminosity (and the mass), and, vice versa, the other 
two models have a fixed value of the luminosity and a var- 
ied value of the temperature. These secondary models are 
used to define the errors on the parameters. In particular, 
we defined the ilcr uncertainty interval as half of the dif- 
ference between the parameters of the quoted "secondary" 
models and the ones of the best fit model. Typically this cor- 
responds to errors of ~ ±12if , ^ zLO.IMq and ~ ±0.01 dex 
in effective temperature, mass and luminosity respectively. 

To analyze the effect of the projection factor p on the 
fit of the radial velocity data, we also tried to vary it in the 
function given by Eq[3] The resulting best fit values of 
the p factor are listed in Tab [2] together with all the other 
structural parameters derived from our analysis. 

We also computed additional models varying the chemi- 
cal composition but these tests did not improve the accuracy 
of the model fitting for none of the investigated Cepheids, 
suggesting the usually assumed Z=0.008 for the metallicity 
of these stars. The results in Tab[2]also show that the value 
of the assumed mixing length parameter (to close the non- 
linear system of dynamical and convective equations) that 
provides the best match is Umi = 1.9 — 2.0 for fundamental 
variables and a slightly smaller value {ami = 1-8) for the 
first overtone one, in agreement with previous theoretical 
results based on the analysis of both Cepheid and RR Lyrae 
prop e rties (see e.g. iBono et al.l I2OO2I: IPi Criscienzo et al.l 



2004 iMarconi fc Clementinil l2005l : iFiorentino et all \2Qm . 
Natale et al. 20081 . and references therein). 

Below, we describe the result of the fit of light and radial 
velocity curves for all the selected Cepheids. 



4.1.1 We 2 

The best fit model selection illustrated above is shown for 
the fundamental pulsator We 2 in Figs lH2l The best fit 
model (central panels in the two figures) has a charac- 
teristic effective temperature of 5925 K and a luminosity 
log(L/LQ) = 3.00 dex. As evident in the top and bottom 
panels, a variation of 25 K in the temperature (at fixed lu- 
minosity) or 0.02 dex in the luminosity (at fixed tempera- 
ture) worsen the match with the data. As mentioned above, 
we define the parameter uncertainty intervals as the half of 
the quoted variations. We also note that the chosen value of 
the projection factor p=1.27 provides an excellent match of 
the model radial velocity with the plotted data, even if the 
best analytic match is obtained for p=1.23. 



4.1 Best fit models 

For each Cepheid we identified a best fit model and other 
four models (hereafter "secondary" models) which are the 
closest (in the sense of the value of the x^ function) to the 
best fit one, according to a reasonablslEl selection of the model 
grid steps (typically 25 K in effective temperature and 0.02 

^ Smaller steps do not produce significant differences in the model 
properties and/or are within the numerical precision of the hy- 
drodynamical code. 



4.1.2 V6 

It is instructive to show the same plots than We 2 for the first 
overtone pulsator V6 (Figs l3l4)l . In this case we note that the 
best fit model is located very close to the first overtone blue 
edge, so that a shift of only 25K towards higher effective 
temperatures (see the top panel of Fig|3} almost quenches 
the pulsation. On the other hand when decreasing the effec- 
tive temperature by 25K (bottom panel of the same figure), 
the predicted amplitudes increase significantly beyond the 
observed ones. If we consider the effect of a variation in the 
stellar mass we note that the corresponding variation in the 
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We2 



—I 1 1 1 — I 1 1 — I 1 1 — I 1 1 — I 1 1 — I 1 — r- 




2.0 



Figure 1. In each one of the three panels on the left, the U, B, V, I and K band observed light curves of the Cepheid We 2 are shown 
from bottom to top (empty squares). The data have been systematically shifted to be shown in the same plot. The radial velocity data 
of the same star is shown in the three panels on the right (empty squares). The solid lines represent the model matched to the data. The 
U band model is plotted with dashed line because it will be excluded in our further analysis (see Sec l4.3l l. The best fit (BF) is plotted in 
the central panels and the models corresponding to — 2ST and + 25T, are shown in the bottom and top panels respectively. 

In the panels showing the radial velocity curves, the models with p-factor free to vary are also plotted (long dashed lines). 



pulsation amplitudes is much smaller, even if non negligible 
when computing the minimization. 

The procedure of best fitting of the radial velocity curve, 
when considering the projection factor as a free parameter, 
provides the value p=1.00, which is significantly different 
from 1.27 and smaller than other typical values adopted in 



the literature. A possible explanation might be a not suf- 
ficiently good quality of the dat a, but we not e that simi- 
lar results have been discussed bv lNatale et~aLl (|2008l ) for S 
Cephei using a well sampled radial velocity curve. 
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Figure 2. The same as Fig. [T]but with models corresponding to log(L/LQ)^ 
2. 



- 2(5(logL/LQ) and log(L/L0) + 251og{L/LQ), for We 



4.1.3 HV 12197 

The two panels of Fig[5] show the best fitting result for 
Cepheid HV 12197. The U band is excluded from our anal- 
ysis due to the poor light curve. According to the best 
fit model this star is the most luminous of our sample 
{\og{L / Lq) — 3.045 dex) and has an effective temperature 
of 5950 K. As for the radial velocity curve, the chosen value 
of the projection factor (1.27) provides a good match with 
the data and it is not significantly different from the value 
p=1.33 obtained from the minimization with variable p. 



4.1.4 HV 12199 

The best fit model for the Cepheid HV 12199 is shown in 
the bottom panels of Fig|6] The luminosity derived from the 
fitting procedure {log{L/LQ) = 2.91 dex) results to be the 
smallest in the selected sample and the effective tempera- 
ture of this star is 6125 K. The model radial velocity curve, 
plotted in the right panel, shows a small discrepancy in its 
amplitude with the data if we fix the value of the p factor 
to 1.27. It is necessary to decrease the projection factor to 
1.17 in order to achieve the best match. It is interesting to 
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Figure 3. The same as Fig. [T]for V6. 



note that, if we focus our attention only on the photometric 
data, the resulting best fit model for IfV 12199 is given in 
the top panels of Fig[Sl with a predicted effective temper- 
ature and log(I//I/0) of 6200 K and 2.93 dex, respectively. 
However, we excluded this model because the radial veloc- 
ity curve, obtained from the chosen value of the projection 
factor, has a too small amplitude to fit the data, (see the 
solid line in the Fig[6]). To account for the amplitude of the 
radial velocity data it would be necessary to decrease the 
projection factor to the too small value p^l.O (dashed line 
in the Figure). This inconsistency that can be indeed due to 
a too high model effective temperature, led us to prefer the 
best model reported in Fig[6] 



4.1.5 HV 12198 

FigEl shows three possible models which describe the pho- 
tometric and radial velocity data of HV 12198. The mini- 
mum for the photometry is obtained for the model at 
Te = 6100 K and log(L/L0) = 3.10 dex (bottom panel). 
However, this does not reproduce accurately the amplitude 
neither of the U band photometry nor of the radial velocity 
data. Decreasing the parameter of the convection efficiency, 
ami , would produce the simultaneous increment of the am- 
plitudes of both light and radial velocity curves. In this way 
we would recover the amplitudes of both the U band and 
the radial velocity curves, but we would decrease the ac- 
curacy of the fit in the remaining bands. As in the case 
of HV 12199, the radial velocity amplitude could be repro- 
duced by decreasing the p factor to the barely acceptable 
value of ~ 1. If we exclude the U band, the model matching 
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Figure 4. The same as Fig. [2] for V6. 



both the photometry and the radial velocity data is the one 
shown in the top panel of FigO corresponding to Te—GOOO 
K and log(L/L0) = 3.1. In this case we obtain an accurate 
fit of the radial velocity data with the chosen projection fac- 
tor value, which decreases to 1.15 if we consider it as a free 
parameter in the minimization. On the basis of these re- 
sults, we decided to consider as our best fit model for HV 
12198 an intermediate case between those already described. 
It is shown in the central panel, includes the U band fitting 
and provides a match of the radial velocity curve which is 
somewhat intermediate between those of the models in the 
top and bottom panels. 



4.2 Comparison with the spectroscopic data 

As a test for the accuracy of our procedure, we compared 
the effective temperature predicted by the best fit models 
for HV 12197 and HV 12199, with the results obtained by 
iMucciarelli et all (|201l[ ) using an independent spectroscopic 
determination. The two panels of Fig|8] show the model 
temperature as a function of the pulsational phase of the 
two quoted Cepheids and their spectroscopic temperatures, 
given by Mucciarelli et al.. (,2011i ) with a typical error bar of 
100 K. 

It is evident that models reproduce the spectroscopic 
data with great accuracy, thus further supporting the pre- 
dictive capabilities of the adopted theoretical scenario, as 
well as of the model fitting technique. 



© 2002 RAS, MNRAS 000, [l]-?? 



8 M. Marconi, R. Molinaro, V. Ripepi, I. Musella and E. Brocato 




Figure 5. The best fit model Ught curves and the model radial velocity curves (solid lines) for the Cepheid HV 12197 are matched writh 
the data (empty squares), in the bottom and top panel respectively. The plotted light curves are in the B, V, I and K bands from bottom 
to top. The best fit radial velocity curve obtained by varying the projection factor is shown as dashed line. 



Table 2. In the top part of the table adopted structural parameters (with uncertainties) of the best fit models obtained for the chemical 
composition Z=0.008 and Y=0.25: (1) Cepheid name, (2-3) mass, (4-5) luminosity, (6) canonical luminosity, (7-8) temperature, (9-10) 
barycentric velocity, (11-12) p factor, (13) convection efficiency parameter. In the bottom part of the table the distance moduli (with 
uncertainties) for all the photometric band obtained from the fit: (1) Cepheid name, (2-3) U band, (4-5) B band, (6-7) V band, (8-9) I 
band, (10-11) K band. 



Structural parameters 

Name ^ <5 (^) log (^) <51og(^) ^(K) ST 7 (Ws) ^7 P Sp a,„, 



HV 12197 4.6 ±0.2 3.045 ±0.012 3.01 5950 ±12 298.3 ±0.9 1.330 Igoog 2.0 

HV 12198 4.2 ±0.1 3.10 ±0.01 2.88 6050 ±12 298.57 ±0.09 1.216 t'o'lH 2.0 

HV 12199 3.5 ±0.1 2.91 ±0.01 2.62 6125 ±12 300.9 ±1.0 1.17 2.0 

We 2 4.30 ±0.15 3.00 ±0.01 2.92 5925 ±12 302.6 ±0.6 1.232 toxnt 1-9 

V 6 4.0 ±0.1 3.03 ±0.01 2.81 6300 ±12 300.6 ±1.0 1.00 : 1.8 



Distance moduli (mag) 

Name fiu Sfiij /ig S/ib /^v ^f^v i^i ^f^i t^K Sfix 



HV 12197 






19.09 


±0.04 


18.96 


±0.05 


18.89 


±0.06 


18.68 


±0.04 


HV 12198 


19.15 


±0.08 


19.13 


±0.05 


18.98 


±0.05 


18.81 


±0.04 


18.60 


±0.03 


HV 12199 


18.97 


±0.05 


18.96 


±0.07 


18.83 


±0.06 


18.68 


±0.04 


18.49 


±0.04 


We 2 


18.70 


±0.08 


18.82 


±0.05 


18.78 


±0.04 


18.71 


±0.03 


18.53 


±0.06 


V6 


19.10 


±0.04 


19.19 


±0.02 


19.03 


±0.02 


18.87 


±0.02 


18.61 


±0.07 
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Figure 6. In the bottom panels the match between the best fit model and data, for the Cepheid HV 12199, is plotted with the same 
meaning of the used symbols than in the Fig[5] The plots in the top panels refer to a model with T=6200 K, log(L/LQ) = 2.93 and 
M/Mq = 3.5. 



4.3 Distance and reddening 

Beyond the intrinsic stellar parameters of the best fit mod- 
els (chemical composition, stellar mass, effective tempera- 
ture, luminosity and convective efficiency parameter) and 
the corresponding projection factor, Tab[2] reports the re- 
sulting distance moduli in all the observed bands, fii (i=U, 
B, V, I, K), with the exception of HV 12197, for which the 
U band data consists only of three measurements and is 
not used to infer the distance. The range of values we find 



for the apparent distance moduli of the selected stars can 
be at least in part understood in terms of differential red- 
dening, but other effects might in principle be at work (see 
below). Using the obatined apparent distance moduli and 
the photometric band e ffective wavelengths, A^, we fitted 
the ICardeUi et al] (Il989l ) extinction law to derive simulta- 
neously the absorption. Ay, and the true distance modulus, 
^0, for the selected Cepheids. To this aim, we minimized the 
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Figure 7. The match between model curves and data is shown for the best fit (central panel) and other two possible models, characterized 
by the same luminosity of the best fit one but slightly lower (Te=6000 K top panels) and higher (Te=6100 K bottom panels) effective 
temperature. As in the previous figures, the photometry is modeled in the left panels and the radial velocity in the right ones. Moreover, 
the light curves in the U, B, V, I and K band arc represented from the bottom to the top. In the radial velocity panels, the dashed lines 
represent the model radial velocity curves obtained from the fit with free p factor. 



following X function by varying the two unknown parame- 
ters: 



fj,i - ( a{xt) + -^^^ ] Av - Ho 



(4) 



whe re the total to selective extinction ratio is fixed to Rv = 
3.3 (|Feast fc Walked 119871 ) . x^ = ^ is the inverse of the i 



band effective wavelen gth and th e expr essions for a{x) and 
b{x) are defined in lCardelU et al.l l|l989l ). 



The results of the fit are shown in fig. |9] and the best fit 
parameters are listed in Tab. O To be conservative, the er- 
rors on the individual band distance moduli include the rms 
of the fitting procedure and the error related to the selec- 
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Figure 9. The result of the best fit of the lCardelli et al.l l ll989t) reddening law is plotted for the five selected Cepheids (solid line). The 
data points (empty squares) represents the distance moduli obtained from the model fitting in the photometric bands considered in this 
work. The two filled triangles represent the U band distance moduli of V 6 and We 2 excluded from the fit. 
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Figure 8. Effective temperature curves predicted by models 
(solid lines) with overp lotted spectroscopic determinations by 
iMucciarelli et alj l l201lf) (empty squares) for HV 12197 and HV 
12199 in the bottom and top panels respectively. 



tion of the best fit modelfl, summed in quadrature. Finally, 
the uncertainties on the parameters of the Cardelli fit have 
been derived from the confidence level at la. In the fit- 
ting procedure, we decided to exclude the U band distance 
moduli of We 2 and V 6 because they significantly deviate 
from the expected trend (see the filled triangles in figH)). 

Inspection of Tab. [3l suggests color excess es- 
timates larger than the typically a dopted value for 
NGC 1866, namely E(B-V)=0.06 mag (|Storm et all l2005l : 
iMolinaro et al.|[2012l . and references therein), with the ex- 
ception of the result for We 2. In fact, the resulting weighted 
mean va lue is E(B-V)=0.11±0.(31 mag , consistent with the 
result bv lGroenewegen fc Salari3 (|2003h . E(B-V)=0.12±0.02 
mag, as based on the simultaneous fit of the NGC 1866 
Cepheid Period-Luminosity relation, in the B, V and I pho- 
tometric bands. 

As for the distance modulus of NGC 1866, a weighted 
mean of the obtained results for the five stars, provides 
fj-o ~ 18.56 ± 0.03 mag. Here, the uncertainty is only sta- 
tistical but we are aware that several systematic effects are 
at work when reproducing observing quantities with pul- 
sation models. First, we have to consider the effect of the 
adopted physical assumptions, namely the equation of state 
and the opacity tables. Previous the oretical investigations 
ijPetroni et al.H2003l : IValle et al.ll2009l ) show that the effect 
of varying these ingredients is marginal and dominated by 
the effect of the model spatial risolution. The latter affects 
the predicted pulsation amplitudes and in turn the intrinsic 
parameters of the obtained best fit models. In particular, in- 
creasing the adopted spatial resolution by 10 mesh zones can 



^ Half the difference between the distance moduli obtained from 
the "secondary" models defined above 



imply a variation of about 100 K in the predicted effective 
temperature and of few hundredths of dex in the predicted 
luminosity level. Another important source of uncertainty 
is the treatment of the pulsation and convection coupling. 
Even if we adopt a nonli near nonlocal time-dependent treat- 
ment of convection (see IStellingwerJI 19821 : iBono et ah I ll999l . 
for details), the equation system is closed by adopting a 
free parameter related to the mixing lenght. Variations of 
the mixing lenght affect the pulsation efficiency and ampli- 
tude. In particular, by varying the mixing lenght parame- 
ters by more than ±0.05 from the value reported in Tab[2j 
we are not able to reproduce the observed curves, with re- 
sulting changes in the predicted luminosity levels smaller 
than ±0.03 dex. Obviously, this is only the effect of varyingh 
the mixing lenght parameter, but within the same turbulent 
convective model. Assuming a different treatment of con- 
vection might in principle produce larger errors even if we 
consider quite encouraging that the application of the model 
fitting technique, from different groups ( e.g. Bono et ajj 
2002; Keller fc Wood 2002; M arconi fc Clc mcntini 20oi; 
[Keller fc Wo'odi I2OO6I ; iMcNamara et all |2007| ) and using 
different approaches to the treatment of the pulsation- 
convection coupling ^ gives consistent resul t s for t he LMC 
distance modulus (|Marconi fc Clementinil l2005l : iMarconil 
i2Q09i). Finally, for what concern the light curves we have 
to consider the uncertainty on the adopted model atmo- 
sphere in transforming bolometric into B,V,I,K variations. 
According to our previous experience we know that theo- 
retical predictions are dependent on the set of static at- 
mosphere models adopted for transforming temperatures 
into colors. For example, changing th e ado pted model at - 
mospheres from ICastelh et all (|1997al lbl) to iKurucj (|l993l ) 
produces colour effects of the order of 0.01 mag. In con- 
clusion, to be conservative, we assume a systematic effect 
of ± 0.1 mag on the inferred distance modulus, as due to 
all the above mentioned theoretical uncertainties. On this 
basis our final estimate of NGC1866 distance modulus is 
/io = 18.56 ± 0.03{stat) ± O.l(syst) mag. 

This result is in agreement within the uncertainty 
interval with the va lue 18.51 ± 0.03 mag, obtained by 
iMolinaro et~aLl (|2012l ) from the Baade-Wesselink method, 
using the same p factor adopted as reference value in this 
work. Their estimates of the distance to HV 12197, We 2 
and HV 12198, are in excellent agreement with the values re- 
ported in Tab. [3] In particular, they found 18.63±0.12 mag, 
18.54±0.09 mag and 18.59±0.08 mag for HV 12197, We 2 
and HV 12198 respectively. For the remaining two stars, HV 
12199 and V 6, they found 18.62±0.10 mag and 18.83±0.11 
mag, respectively, both systematically longer than our esti- 
mates, although consistent with them within the errors. In 
a recent wor k , usin g the infrared surface brightness method, 
IStorm et al. I (|201ll ) obtained the distance for a sample of 
Cepheids in the LMC, including HV 12197, HV 12199 and 
HV12198, and their final value (18.45 ± 0.04) is in agree- 
ment with our results within the errors. However, they 
used a period dependent p factor given by the equation 
p = 1.550(±0.04) - 0.186(±0.06) logP. This is required to 
obtain distances to LMC Cepheids independent of their pul- 
sation periods and distances to Galactic Cepheids in agree- 
ment with th e HST para llaxes, with their surface brightness 
method (see IStorm et al... 2011). As stated by the authors 
themeselves, this relation is not easily reconciled with recent 
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Table 3. Parameters obtained by fitting Cardelli extinction law: 
the Cepheids are listed in the first column, the second and third 
columns contain, respectively, the extinction, Av , and the red- 
dening, E(B — V), the intrinsic distance modulus, is in the 
last column. 

Name Ay (mag) E(B-V) (mag) /xq (mag) 



HV 12197 0.33 ±0.07 
HV 12198 0.43 ± 0.07 
HV 12199 0.37 ±0.07 
We 2 0.23 ±0.10 
V 6 0.46 ± 0.07 



0.10 ±0.02 
0.13 ±0.02 
0.11 ±0.02 
0.07 ±0.03 
0.14 ±0.02 



18.65 ±0.06 
18.55 ±0.05 
18.46 ±0.06 
18.55 ±0.08 
18.57 ±0.06 



theoretical work (e.g. Nardetto et al. 2009) and provides p- 
factor values for short-period Cepheids (not less than 1.4), 
significantly larger than the results derived in the present 
study. This discrepancy could be due, at least in part, to 
limitations of our treatment of the coupling between pul- 
sation and convection (see discussion above) but we have 
to note that the debate on the p factor and on its possi- 
ble dependence on the p ulsation period is s till open in the 
recent literature (see e.g. iNgeow et al.ll201^ . and references 
therein) . 



4.4 Mass— Luminosity relation 

The obtained stellar masses for the Cepheids in NGC 1866, 
as reported in Tab[2] cover a range of values that might be 
the signature of differential mass loss . This occurrence is i n 
agreement with previous findings by iBrocato et al. I (|2004l ). 
Finally, we compared the derived masses and luminosities, 
as reported in T ab[2l with an evolut ionary Mass-Luminosity 
relation (MLR) jBono et al.ll2000bl ). either neglecting or in- 
cluding mild oversh ooting according to the prescriptions by 
IChiosi et all (Il993l ). 

Fig llOl clearly shows that the analyzed Cepheids do not 
follow one of the two relations, but they are randomly placed 
at intermediate luminosities between those predicted by the 
canonical and the mild overshooting MLRs, with the ex- 
ception of the HV 12199 which results to be slightly more 
luminous than the mild overshooting prediction, although 
consistent with it within the uncertainties. 

The fact that different luminosities are predicted for 
a given mass might suggest that the investigated Cepheid 
do not follow a MLR but are instead stochastically affected 
by some noncanonical phenomenon, likely a combination of 
mild overshooting and mass loss. However the application of 
the method to a larger sample of pulsators is needed in order 
to draw any reliable conclusion on the cause of the overlu- 
minosity distribution with respect to the canonical one. 



5 CONCLUSIONS 

We have used nonlinear convective pulsation models com- 
puted by our team to reproduce the multifilter (U,B,V, I and 
K) light and radial velocity curves of five Classical Cepheids 
in NGC 1866, a young massive cluster of the Large Mag- 
ellanic Cloud. The resulting best fit models give us infor- 
mation on the intrinsic stellar parameters and the individ- 




0.55 0.60 0.65 

log(M/Mj,) 



Figure 10. The masses and luminosities of the analyzed 
Cepheids (empty squares) are compared with the canonical Mass- 
Luminosity relation (solid line) and with the mild overshooting 
Mass— Luminosity relation (dashed line). 



ual distances of the investigated Cepheids. In the case of 
HV 12197 and HV 12199 the obtained effective tempera- 
ture and its variation with the pulsation phase has been 
found to be in very good agreement with the spectroscopic 
determinations within their uncertainties. The masses and 
luminosities, obtained for all the five investigated pulsators, 
from this model fitting technique satisfy a slightly brighter 
Mass-Luminosity relation than the canonical evolutionary 
one, indicating that noncanonical phenomena such as mild 
overshooting and/or mass loss are at work. As for the in- 
ferred distances, the individual values have been found to be 
consistent with each other within the uncertainties. More- 
over, their weighted mean value corresponds to a distance 
modulus of 18.56±0.03 (stat) ± 0.1 (syst) mag, in agree- 
ment with several independent results in the literature. In 
particular the obtained result for the distance to NGC1866 
is in excellent agreement with the results obtained from 
the application of the model fitting techni que to LMC field 
Cepheids, RR Lyrae and & Scuti variables (iBono et al.|[2002l : 
iMarconi fc Clementinill2005l : iMcNamara et al.ll2007l ). 
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